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Abstract 

We numerically solve fully (3+l)-dimensional relativistic hydro-dynamical equation coupled with the 
baryon number conservation law without spatial symmetry. We discuss the effect of transverse expansion 
based on the deviation our numerical result from Bjorken's scaling solution. We analyze the space-time 
evolution of the QGP gas in the case of non cylindrical initial conditions. 

1 Introduction 

The various kinds of collective flow phenomena such as directed flow, elliptic flow and radial flow has 
been observed in recent experiments at AGS [El and SPS 0. It is a matter of interest that such flow are 
results of hydrodynamical motion of hadronic fluid. Our first trial to tackle the problem is to develop 
(3+l)-dimensional hydrodynamical model. Assuming the local thermal equilibrium for hot and dense fire 
ball produced in ultra relativistic nuclear collisions, we analyze the evolution of the fire ball based on the 
(3+l)-dimensional hydrodynamical model. The hydrodynamical model for Quark-Gluon Plasma (QGP) 
fluid has already been discussed in many papers since Bjorken first introduced the simple scaling model 
based on (l+l)-dimensional expansion picture ||. For simplicity, cylindrical symmetry is assumed in 
usual hydrodynamical ^, [||^] analysis, but this assumption disables us from discussing the anisotropic 
collective flow. In this paper, in order to investigate collective flow not only in the central collisions but 
also in the non-central ones we numerically solve the (3+l)-dimensional relativistic hydrodynamical 
equation coupled with the baryon number conservation law. 

2 The relativistic hydrodynamical model 

The relativistic hydrodynamical equation for perfect fluid is given as 

9 M T^ = 0, (1) 

where T M " is energy momentum tensor, 

T^ v = ell^U" — P(ff M " — U^U"). (2) 

Here, e is energy density, P is pressure, metric tensor is g^" = diag.(l, —1, — 1, —1) and local velocity is 
£/ M = (l,v x ,Vy,v z )/~/ respectively. In order to take account of the finite baryon number density, we must 
solve the baryon number conservation law, 

du{n B {T, l x)U"} = Q, (3) 

also, where ns(T, /i) being baryon number density. Through the time like component of Eq.(Q), Uud^T^" = 
0, Eq.(^) and thermodynamical relation, e + P = TS+fins, we can obtain the conservation law of entropy 
density current, S" = S£/ M , 

= 0. (4) 
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Our numerical algorithm solving the hydrodynamical equation is based on the entropy conservation law 
Eq.(§. 

In order to solve the hydrodynamical equation, the equation of state is needed. Though we consider 
the QGP gas and the hadron gas(excluded volume model ||) for the realistic model equation of states, 
in this paper we adopt the QGP gas for numerical simplicity. The QGP gas model of massless Nf flavor 
quarks is given by, 

32 + 21iVf 2 4 N f a a N f 4 
P = 18Q * T + —» q T + — N , (5) 

where the number of flavor being Nf = 3 and chemical potential for quarks being fi q = /i/3. 

3 The numerical calculation 

We solve the (3+l)-dimensional hydrodynamical equation without symmetrical conditions by using an 
algorithm in which lattice points of volume element is moved along local velocity and the entropy con- 
servation law Eq.(^) is adopted explicitly. D. H. Rischke et al. discuss relativistic hydrodynamics in 
(3+l)-dimensional situation and collective behavior by using Eulerian hydrodynamics Q. Our numeri- 
cal calculation is explained briefly as follows: 

In the first step the coordinates x m = X m (t, i, j, k) (m — 1, 2, 3) of lattice points at time t + At are 
replaced by 

X m (t + At, i,j, k) = X™ (t,ij, k) + ^ (6) 

In the determination of lattice points in Eq.(^|), the coordinates move in parallel with nsll^, SIP 1 . 
In the next step the local velocity is determined by, 

v m (t + At,i,j,k) = v m (t,i,j,k)+d t v t (i,j,k,t)At (7) 

3 

+ (*. 3, k, t) {X m (t + At,i,j,k)~X m (t,i, j, fc)) 

where 9 M u M obtained from Eq.(jl|), Eq.(|^) is used. 

In the final step the temperature and chemical potential of lattice points is calculated by using Eq.(^|), 
Eq.(§). 

4 Comparison with Bjorken's solution 

Comparing our numerical solutions with Bjorken's scaling solution v z — z/t, we can easily evaluate the 
effect of the transverse flow. Based on Bjorken's scaling solution and Eq.M) entropy density is given as, 

S(t) = S(toA (8) 

T 



where proper time t, t = \/t 2 — z 2 . In order to make comparison clear, in this section the velocity 
of our model in the longitudinal direction is fixed to the Bjorken's scaling solution. We put the initial 
temperature distribution and chemical potential distribution respectively as follows: 



fi(t ,x,y,z) = fj, < exp 



(z — Zb) 



2a B 



+ exp 



(z + z B ) 2 
2a b 2 



{-feS^-«>-w e <»-»>}' (I0) 

where x — y = z = 1.0 fm, z B = 0.7 fm, a x — a y = <j z — 1.0 fm, a B = 0.7 fm, T = 200 MeV, 
(j,q — 210 MeV and the initial transverse velocity is set to 0. We focus to the volume elements at 
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Figure 1: the comparison with Bjorken's 
scaling solution 




Figure 2: the difference between our so- 
lution and Bjorken's scaling solution by 
changing xq, yo from 1.0 fm to 6.0 fm 
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Figure 3: the non cylindrical initial condi- 
tion of temperature distribution at z = 
fm 
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Figure 4: The solid line stands for freeze- 
out which we assume from chemical freeze- 
out (the dashed line) and thermal freeze-out 
(the dotted line). 



(x,y, z) — (0,0,0) for the comparison of our numerical solution with Bjorken's scaling solution. Figure 
1 shows that our numerical calculation is coincident with Bjorken's scaling solution up to r = 3.0 fm. 
After r = 3.0 fm the difference between them increases with the proper time because of the transverse 
flow. We define a characteristic time tb at which instance the difference between two models becomes 
larger than 0.1 % . Figure 2 indicates that tb is almost proportional to the initial xo, yo during 1.0 fm 
to 6.0 fm, and Bjorken's scaling solution seems to be a proper solution for a large system. 



5 Non cylindrical initial conditions 

In this section we investigate the space-time evolution of the hydrodynamical flow with non cylindrical 
initial temperature distribution and chemical distribution. Other parameters are put as the previous 
section. Figure 3 shows the initial temperature distribution in b = 0.8 fm and 72 = 1.0 fm. We evaluate 
the effect of the flow in particle distributions by giving initial conditions like this, though this condition 
is not realistic to analyze the experimental data. We use Cooper-Frye formula fl0| for particle emission 
from hadronic fluid, 

dN _ g h f M 1 

VP ~ (2tt)3 J a M exp[(P„[/" - n)/T[] ± 1 (LL> 

for evaluating one-particle distributions. We assume that hadronization process occurs when the tem- 
perature and chemical potential in the volume elements cross the boundary (the solid line) in fig.4. The 
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Figure 5: the space-time evolution of the flow at z = fm (b = 0.8 fm) 



1.1 



o 1.0 



0.9 







---- b=0 [fm] b=0.4 [fm] — 


b=0.8 [fm] 



> 10 

q 1 

t> 10 
%10 



H£io 




100 200 300 

azimuthal angle (<]> ) cleg. 

Figure 6: the azimuthal distribution in 
changing b 
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Figure 7: the dependence of the transverse 
momentum distribution on the flow (b = 0.8 
fm) 



solid line in fig. 4 is so designed that the freeze-out temperature becomes 140 MeV at vanishing chemical 
potential, based on chemical freeze-out and thermal freeze-out model in ref. ^^|. Several calculations 
are made for different initial conditions. Figure 5 displays the space-time evolution of the flow in b — 0.8 
fm. Figures 6 and 7 show the azimuthal fluctuation of particle number and Pt distribution which are 
caused by non-cylindrical properties of transverse expansion. Figure 6 indicates that the variation in the 
azimuthal distribution increase as separation of two initial blobs increases. Figure 7 indicates the influ- 
ence of the flow increase with Pt- The yield at <j) = 90°, 270° is large, because freeze-out hypersurface is 
large in these directions as fig. 5 shows. Furthermore the transverse momentum distribution at <j) = 90°, 
270° is flatter, since the flow is pushed out at = 90° , 270° . 



6 Summary 

We solved (3+l)-dimensional relativistic hydrodynamical equation without cylindrical symmetry condi- 
tions by Lagrangian hydrodynamics. We discussed the effect of the transverse flow and confirmed that 
Bjorken's scaling solution is a proper solution in a large system by making a comparison with numerical 
calculation. The effect of the flow to the particle distributions was also investigated. The influence of 
flow is large at <j> = 90°, 270°. This is because the freeze-out hypersurface is large at <j> = 90°, 270° and 
the flow is pushed out in the direction of <j> = 90° , 270° . 

We need to use more realistic initial conditions in the analysis of the experimental data. We plan 
to adopt the output from the event generator as initial conditions and to use the equation of states 
including phase transition from the QGP phase to the hadron phase. Investigating the collective flow in 
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experimental data is our next task. 



References 

[1] J. P. Wessels: E877, Nucl. Phys. A638(1998), 69c; H. Liu : E895, Nucl. Phys. A638(1998), 451c; 
S. A. Voloshin: E877, Nucl. Phys. A638(1998), 455c 

[2] H. Appelshauser et al. : NA49, Phys. Rev. Lett. 80(1998), 4136; M. M. Aggarwal et al.: WA98, 
|nucle-ex/9807004 

[3] J. D. Bjorken, Phys. Rev. D27 (1983), 140. 

[4] T. Ishii and S. Muroya, Phys. Rev. D46(1992), 5156 

[5] J. Sollfrank, P. Huovinen, M. Kataja, P. V. Ruuskanen, M. Prakash, R. Venugopalan, Phys. Rev. 
C55(1997), 392 

[6] C. M. Hung and E. Shuryak, Phys. Rev. C57(1998), 1891 

[7] B. R. Schlei and D. Strottman, Phys. Rev. C59(1999), R9 

[8] D. H. Rischke, S. Bernard, J .A .Maruhn, Nucl. Phys. A595(1995), 346; D. H. Rischke, Nucl. Phys. 
A610(1996), 88c 

[9] D. H. Rischke, M. I. Gorenstein, H. Stocker, and W. Greiner, Z. Phys. C51(1991), 485 

[10] F. Cooper and G. Frye, Phys. Rev. D10(1974), 186 

[11] U. Heinz, Nucl. Phys. A638(1998), 357c 



5 



